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Abstract 

Energy efficiency has become an encouragement, and more than this, a requisite for the design of 
next-generation wireless communications standards. In current work, a dual-hop cognitive (secondary) 
relaying system is considered, incorporating multiple amplify-and-forward relays, a rather cost-effective 
solution. First, the secondary relays sense the wireless channel, scanning for a primary network activity, 
and then convey their reports to a secondary base station (SBS). Afterwards, the SBS, based on these 
reports and its own estimation, decides cooperatively the presence of primary transmission or not. In 
the former scenario, all the secondary nodes start to harvest energy from the transmission of primary 
node(s). In the latter scenario, the system initiates secondary communication via a best relay selection 
policy. Performance evaluation of this system is thoroughly investigated, by assuming realistic channel 
conditions, i.e., non-identical link-distances, Rayleigh fading, and outdated channel estimation. The 
detection and outage probabilities as well as the average harvested energy are derived as new closed- 
form expressions. In addition, an energy efficiency optimization problem is analytically formulated and 
solved, while a necessary condition in terms of power consumption minimization for each secondary 
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node is presented. From a green communications standpoint, it turns out that energy harvesting greatly 
enhances the resources of secondary nodes, especially when primary activity is densely present. 

Index Terms 

Cognitive relaying systems, cooperative spectrum sensing, detection probability, energy efficiency, 
energy harvesting, green communications. 

I. Introduction 

Recently, cognitive radio (CR) has emerged as one of the most promising technologies to 
resolve the issue of spectrum scarcity, caused by the escalating growth in wireless data traffic of 
next-generation networks [Q], [j2|. One of the principal requirements of CR is the effectiveness 
of spectrum sharing performed by secondary (unlicensed) nodes, which is expected to intelli¬ 
gently mitigate any harmful interference caused to the primary (licensed) network nodes. This 
requirement is directly related to the accuracy of spectrum sensing/sharing techniques, reflecting 
the reliable detection of primary transmission. Moreover, to further guarantee a sufficient quality 
level of primary communication, the transmission power of CR is generally limited, such that 
its interference onto primary users remains below prescribed tolerable levels. However, this 
dictated constraint dramatically affects the coverage and/or capacity of the secondary communi¬ 
cation. Such a condition can be effectively counteracted with the assistance of wireless relaying 
transmission. In particular, the rather feasible dual-hop multi-relay communication scheme with 
best relay selection is of paramount interest due to its enhanced performance gains (e.g., see 
[[3J-[l9l and references therein). 

Building on the aforementioned system deployment, the spectrum sensing process can also 
be significantly enhanced. By means of the so-called cooperative sensing (e.g., see IflOl . ifTTTO . 
each secondary node may sense the channel in fixed time periods and then forward its sensing 
measurement to a central secondary base station (SBS), which acts as a fusion center. The latter 
entity is responsible for the final decision on a primary transmission occurrence, benefiting from 
the spatial diversity of several sensing reports. Such a distributed (cooperative) spectrum sensing 
was shown to deliver much more accurate decisions than local (standalone) sensing regarding 
the detection of primary transmissions iflTTl . 

On another front, driven by the ever increasing economical and environmental (e.g., carbon 
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footprint) costs associated with the operating expenditure of communication networks, energy 
efficiency (EE) has become an important design consideration in current and forthcoming wireless 
cognitive infrastructures lfT2ll - IIT4l . Taking into account that nodes in wireless communication 
systems are mostly battery-driven while governed by power-hungry equipment, EE appears to 
be an essential requirement. 

A. Related Work and Motivation 

In cooperative CR systems, two different types of relaying protocols have dominated so far, 
namely, amplify-and-forward (AF) and decode-and-forward (DF) Ifl5l . It is noteworthy that AF 
outperforms the computational-demanding DF in terms of EE and/or power savings. In ffl6l 
and ifTTIl . the tradeoff between consumed power and performance of the secondary system was 
investigated. However, in these works, the objective was performance enhancement (in terms 
of either outage probability |[76ll or throughput ifTTlO not EE and/or minimization of power 
consumption. In Ifl8l - lf20l . the authors focused on the energy minimization of cognitive relaying 
networks. Yet, EE and data transmission requirements (e.g., in terms of a data rate and/or error 
rate target) were not jointly considered in these works. In lfl3l . the latter problem was jointly 
considered, given a predetermined detection probability on the primary nodes activity (i.e., 
providing conditional expressions). A recent approach dealing efficiently with power savings 
and minimization of consumed energy is energy harvesting ETI . Powering mobile devices by 
harvested energy from ambient sources and/or external transmission activities renders wireless 
networks not only environmentally friendly but also self-sustaining. From the CR perspective, 
lf22l . Il23ll investigated cognitive systems enabled with energy harvesting equipment, but coop¬ 
erative spectrum sensing or relayed transmission was not considered. 

In this paper, we focus on the EE optimization of CR cooperative systems, while an op¬ 
portunistic strategy incorporating energy harvesting is introduced. According to the proposed 
strategy, secondary nodes switch between data transmission and energy harvesting depending on 
their sensing decision on the existence of primary nodes activity. More specifically, a dual-hop 
relaying system with multiple relays is adopted for the secondary system, where the end-to- 
end (e2e) communication is facilitated via a best relay selection policy. AF relaying protocol is 
implemented via semi-blind relays and fixed gains, a rather cost-efficient solution. Cooperative 
spectrum sensing is performed in a fixed sensing time duration, followed by a reporting of relay 
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sensing measurements to SBS. Upon the aggregate decision at SBS, secondary nodes enter into 
either the energy harvesting phase (if primary transmission is detected) or the transmission phase 
(if primary transmission is not detected). In order to avoid unexpected co-channel interference to 
primary nodes (e.g., during reporting, transmitting and backhaul signaling from SBS to relays), 
the average interference temperature constraint is considered. 

B. Contributions 

The main contributions of this paper are summarized as follows: 

• A novel hybrid mode of operation is introduced for the CR cooperative system, where the 
secondary nodes opportunistically switch between data transmission and energy harvesting, 
upon the detection of primary activity. 

• A new exact closed-form expression for the detection probability of the considered con¬ 
figuration is presented, assuming that all the involved channels undergo independent and 
non-identically distributed (i.n.i.d.) Rayleigh fading conditions (i.e., different link distances 
among the nodes, appropriate for practical applications). 

• A new exact formula for an important system performance metric is derived, namely, the 
outage probability, when the secondary system enters into the data transmission phase. 
The practical scenario of outdated channel state information (CSI), between the reporting 
and transmission phases, is also considered. Moreover, an exact formula for the average 
harvested power is presented, when the secondary system enters into the corresponding 
phase. 

. Based on the aforementioned performance analysis, the overall energy consumption at each 
secondary node is formulated. An optimization problem aiming at minimizing this energy is 
introduced, whereas a necessary and sufficient optimality condition is manifested. Finally, 
based on the enclosed analysis, the optimal sensing time is numerically evaluated given the 
required constraints. 

C. Organization of the Paper 

In Section [TT] the signal model and the proposed opportunistic strategy for energy efficient cog¬ 
nitive relaying systems is presented. Subsequently, closed-form expressions for key performance 
measures of the secondary system are included in Section Hill Based on these expressions, 
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Fig. 1. The considered system configuration with a secondary source (S), M relay (R) nodes and an SBS, i.e., a destination 
node (D), under the presence of L primary (P) nodes. 


the energy consumed from each secondary node is analyzed in Section |IVl followed by the 
introduction of an optimization problem on minimizing the overall energy consumption of the 
CR cooperative system. Indicative numerical results and useful discussions are presented in 
Section [V] while some concluding remarks are presented in Section EO 

D. Notation 

Throughout this paper, the following notations are used: E[-] stands for the expectation operator 
and Pr[-] returns probability. Also, fx(r) and F x (-) denote, respectively, probability density 
function (PDF) and cumulative distribution function (CDF) of random variable (RV) X. Fur¬ 
thermore, r(-,-) represents the upper incomplete Gamma function ||24l Eq. (8.350.2)], Ei(-) is 
the exponential integral li24l Eq. (8.211.1)], J 0 (-) denotes the zeroth order Bessel function of the 
first ki nd ll24l Eq. (8.411)], 7 0 (-) is the zeroth order modified Bessel function of the first kind 
|24l Eq. (8.431)], and K n (•) is the nth order modified Bessel function of the second kind lt24l 
Eq. (8.446)]. 


II. System Model 

A secondary (cognitive) dual-hop system consisted of a source (S') communicating with a 
destination ( D) via M relay (/?, with 1 < i < M) nodes is considered!] Direct communication 
between source and destination is not available due to the long distance and strong propagation 

’Note that the terms D and SBS will be interchangeably used in the rest of this paper. 
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Fig. 2. Frame structure of the proposed approach. 


attenuation, while keeping in mind that in secondary systems the transmission power must in 
principle be maintained in quite low levels. The system operates in the vicinity of another 
licensed primary network, which consists of L primary (Pj with 1 < j < L) nodes (c.f., Fig. [Cl¬ 
in current study, we assume that all the involved signals are subject to i.n.i.d. Rayleigh fading 
as well as additive white Gaussian noise (AWGN) with a common power N 0 . Thus, PDF of the 
instantaneous signal-to-noise ratio (SNR) is given by 


UA X ) 


N 0 exp 


/ Ngx \ 


Pij'yij 


x > 0 


( 1 ) 


where 7 ^, p i} j and 7 it j = ddj 1 ’ 3 denote the instantaneous SNR, the signal power and the average 
received channel gain, respectively, from the zth to jth node. Moreover, d l , J and a hl represent 
the corresponding distance and path loss factor, respectively. Usually, a hJ e {2, 6 } denotes 
free-space loss to dense urban environmental conditions, correspondingly. 


A. Protocol Description 

The basic lines of the proposed approach are sketched in Fig. [3 The secondary nodes operate 
in a time division multiple access scheme, where sensing and transmission or harvesting phases 
are periodically alternating. 

1) Sensing Phase: First, the relays and the destination enter into the sensing phase where they 
listen to the presence of primary users’ signals over the shared spectrum band within a fixed 
sensing duration T s . We assume that the sensing time is smaller than the channel coherence 
time of primary transmission, such that the magnitude of the channel response remains constant 
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within this phase. This is a reasonable assumption since the sensing duration is, in principle, 
maintained quite small (e.g., see lf25l ll26ll and references therein). The received signal at the 
/th relay or destination can be expressed as 

L 

yp,x = ^2 e jy /p-pg Pj , xSj + n x , 1 < i < M, X G {R t , D} (2) 

3 = 1 

where 9j = 1 or 0 when the jth primary signal is present or absent, respectively. Also, p p , fjp rX , 
Sj and n x denote the transmit power of primary nodesthe instantaneous channel gain from Pj 
to X, the transmitted data of the jth primary node and AWGN at X, respectively. In particular, 
9j follows the Bernoulli probability mass function (PMF), such that = 1 or 9j — 0 with 
probabilities Pr [9j = 1] and (1 — Pr \9j = 1]), respectively. 

2) Reporting Phase: Next, each relay amplifies and forwards its local sensing measurement 
to the destination on its particular time slot (which is a priori reserved from the system), entering 
into the reporting phase. Hence, the received signal at the destination, forwarded by the Ah relay 
at its allocated time slot, yields as 

VRi,D = JpZWdGk,*™ + n D (3) 

where P'^o an d 9n, d are the transmission power and channel gain from R, to D, respectively, 
whereas G-R,i denotes the fixed gain of the Ah relay, all indicating the reporting phase. 

In the sensing phase, each relay performs spectrum sensing on a per sample basis and then 
forwards its report to SBS. Thereby, defining that U samples are provided within the sensing 
phase duration, the following procedure occurs: a) secondary relays sense the spectrum for 
sample-time 1, then go to reporting phase for sample 1; b) they sense the spectrum for sample¬ 
time 2, then go to reporting phase for sample 2; and so on for U samples. Notice that the 
latter alternation between sensing and reporting phases is a natural outcome of the AF relaying 
protocol’s structure (receive and forward without further processing), whereas it further enhances 
the joint spatial-time diversity for the spectrum sensing process (since U > 1). 

3) Harvesting Phase: Then, the destination determines the presence of a primary transmission 
or not, according to the received signals’ power. In fact, it compares the maximum from M + 1 
signals (from the relays and its own) with a predetermined power threshold value A. In the case 


2 Without loss of generality and for the sake of clarity, a common power profile for the primary nodes is adopted. 
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when this signal power is greater than A, a detection event is declared in a subsequent time 
slot and all the relays initiate a harvesting phase, collecting energy from the occurring primary 
transmission(s). To this end, we have that 


Ph = Pr 


max < 7 p,D) max 


(iz) 

%2 Li 


M 


|=1 


> A 


u=1 


( 4 ) 


where Pd stands for the detection probability, while represents the e2e SNR at the sensing 
phase from the (potentially active) primary nodes to destination via the ith relay. For analytical 
tractability, we assume independence amongst the samples, which, nonetheless, represents a 
common yet efficient approach (e.g., see lt28l . If29ll and references therein). Moreover, the 
harvested power at the /'th relay conditioned on a detection event, reads as Il30l Eq. (2)] 

L 

E H ,i = vPdPp ^2 QpPj,Ri 1 9Pj,Ri | ( 5 ) 

3 = 1 


where p e (0,1] is the radio frequency-to-direct current (RF-to-DC) conversion efficiency. 

4) Transmission Phase: On the other hand, if the power of the signal in © is lower than A, 
primary transmission is not detected with a probability 1 — P d . In such a case, capitalizing on the 
status of channel gains from all the relay-to-destination links (collected from the aforementioned 
reporting phase), the destination selects the relay with the highest instantaneous SNR (i.e., R s 
with s determined by the condition 7 p s d = maxi{'y RhD }f£ 1 ) and broadcasts this information in 
the subsequent time slot. In turn, the selected relay informs the secondary source to enter into 


the transmission mode of operation^ Based on this call, the secondary system enters into the 
transmission phase, while the source communicates with the destination via the selected relay 
(all the other ones stay idleQ). 

In fact, the classical half-duplex dual-hop AF relaying protocol is established at this stage, 
where the source-to-relay and relay-to-destination links occur in orthogonal transmission phases 
(e.g., in two consecutive time slots). Hence, the received signals at the relay and destination are, 


3 The event of no signal returning from any relay back to the source, at this stage, can be interpreted as a triggering of 
harvesting phase for the secondary source. 

4 In current study, we focus on single-band (low-complex/low-cost) AF relay nodes. In addition, we assume that all active 
transmissions in the vicinity of the secondary system are realized by either the primary or the secondary nodes in that frequency 
band. 






respectively, given by 


r S,R, 



hs,R s z + n R 


( 6 ) 


and 



( 7 ) 


f'T) 

where z, hs, Rs , P RsD and G r .s correspond to the source data, instantaneous channel gain 


from S to R s , transmission power from S to R s , and fixed gain of the sth relay during the 


transmission phase, respectively. Also, h Rs:D denotes the channel estimate of the selected relay 


to the destination, based on the instantaneous channel status derived from the previous reporting 
phase. It is noteworthy that h RsyD could vary from the actual h Rs)D due to a possible outdated 
CSI at the destination. This condition is realized when a feedback delay and/or rapidly varying 
fading channels between the reporting and transmission phases are present. As such, the channel 
estimate is formed as m 



( 8 ) 


where w Rs .n is a circularly symmetric complex Gaussian RV with the same variance as h Rs .n, 
while p s denotes the time correlation coefficient between h RstR > and h Rst o, defined as ll32l Eq. 
(4.1-63)] 



where fP and are the maximum Doppler frequency on the i—j link and the time difference 


between the actual instantaneous channel status and its corresponding estimate, respectively. It 
follows that in the case when the channel instances remain constant between the reporting and 


transmission phase, then p = 1 and h Rs}D = h RgyD . 


B. Transmission Power of Secondary Nodes 

Although the transmission power of the primary service takes arbitrary values, this condition 
does not apply for the secondary system. We adopt an average interference constraint for the 
transmission power of secondary nodes, taking also into consideration the maximum output 
power, namely, P max . Thereby, for the reporting phase, the following condition should be satisfied 



(10) 
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where Q represents the received power threshold that should not be exceeded at the primary 
nodes, and 


q Ri = max{ 7 iiP y =1 . 


In the case when the system enters into the transmission phase, we have from 
that 

Q 


and 


p S ,R.=m m|P_ Tl _ p<i)E[(B] 

( T) ■ f D Q 

P "-" = ,,1, T“FWW 

where q s = max,- {7 s.Pj}^ and q Rs = max{ 7 Rs:P) }J J =1 . 


(ID 

Eq. ( 6 )] 

( 12 ) 

(13) 


III. Performance Analysis 

We commence by describing the instantaneous e2e received SNR at the destination in the 
aforementioned phases of the proposed strategy. Then, the detection probability is derived, 
followed by some important system performance measures for the secondary system, namely: i) 
the average harvested power of each relay during the harvesting phase, and ii) the e2e outage 
probability during the transmission phase. 


A. SNR Statistics in the Sensing Phase 

Based on d2]). the received SNR from primary nodes to D involves the sum of i.n.i.d. expo¬ 
nential RVs (with different link distances) and is obtained bjf] 

7 p,d--tt^ 2 OjlPj, d- ( 14 ) 

N o j=1 

According to the total probability theorem, PDF of 7 P D is derived as 

L 

f^A x ) = 

r =1 



= r 


flP,D\r{ x ) 1 


(15) 


5 Since we model the signals as RVs with known transmission power, energy detector is adopted at the receiver, as being the 
optimal technique to detect the primary transmission(s) OH 
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where Pr[(]TV[ 0 j = 1 ]) = r] denotes the probability that r primary nodes are active (transmitting), 
given by 


Pr 


= i] 


= r 


fr= (^Pr[% = ini-Pr[», = l]) L -', 


while f lPD \ r {x) denotes the conditional PDF of 7 pp given r primary signal transmissions, which 
reads as If34l Eq. (5)] 


flP,n\r( X ) = l Y. 


n 

j —1 


\ 


lP k ,D ~ lPj,D 


exp 


NqX 


PplP k ,D 


Therefore, the unconditional PDF of 7 PD is expressed as 


/ 


Nn 


■ "P ' r =1 k =1 


\ 


n 


?= 1 7 P k ,D ~ lPj,D 


exp 


N 0 x 


PplP k ,D 


Moreover, CDF of 7 ^ is directly obtained as 

( 


f ipA x ) = 1 - 


r=l fc=l 


\ 


n 


= 1 7 P k ,D - lPj,D 


\j¥=k 


1Ph,d exp 


/ 


iV 0 a: 


PplP k ,D 


B. SNR Statistics in the Reporting Phase 

Based on (O, the SNR of P — Ri — D link (with 1 < i < M ) is given by 

„, = 7 !?^ 

7e2e '‘ 


where 


7i 


(TV) _A_ Pp 




(TV) 

(TV) a P RiP„,(K) 


No — -• ■ iv 0 

j=i 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


The parameter = 1/ jiVo) indicates a constant parameter, which is related to the value 
of fixed gain of the ith relay for the reporting phase. Among some popular precoding designs 
for this parameter, there is quite an efficient one |35|. yielding 

1 


U (n) = I E 


7!? +1 


( 21 ) 
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Lemma 1: A closed-form expression for the CDF of 7 ^ ?: under i.n.i.d. Rayleigh fading channels 
reads as 

( , \ 


with 


7<k>) = i - £/-£ 


2 Af» fPPPPP 


x exp 


r=l y/Pi(p ( IU,D'yRi,D)* 

N 0 x 


n 


j =i ,jfe 

V?¥fc 


7p fc ,r> — 7 p,,d 


PplP k ,Ri 


Ki 2N (U 


U™x 


\ (K) - 

\ PpPRi,DlPk,RilRi,D 


«7= £/,£ 


r=l fc=l 


c 

n 

j=l 

Vv* 


\ 


lP k ,D ~ lPj,D 


r (o, 


N 0 


PplP k ,Ri 


A" 


/ 


exp 


JVo 


Pp7P k ,Ri 


( 22 ) 


(23) 


and 


(U) 

PRi,D ~ 


P 


1 gfajj 

0 


-1 


L L L 


(-1 )* 


e^hXXX-'-X- 7 . . v 

l= i fc=o " i=1 v + ^ =i 

Tl 1 ^ • ^71 ^ ^ fc 

Proof: The proof is relegated in Appendix |Al 


(24) 


(25) 


C. Detection Probability 

Capitalizing on the aforementioned statistics, P d can be obtained in a closed-form expression. 
Proposition 1: Detection probability for a given power threshold A, P d (A), is given by 


Pd( A) = 1 - 


M 




i =1 


U 


= l-A(Af, 


(26) 


where 
A(A) = 


M 

n 

i =1 


/ 


X/ Y, 2No ^^ Pk,RM ' n)x 
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y^iP&'TRi.D^ 


r =1 
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J= i 7 p*,p> - 7Pj,o 


Ki 2N 0 


u^h 

(1Z) - 

ppPR i ,D'yp k ,R i 'yR i ,D 
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exp 


( Nq\ 
\PpiP k ,Ri 
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X 


1 -'£fr'E 




r =1 


k =1 


\ 


n 


= 1 7 Pk,D - lPj,D 


\j¥=k 


7 P k ,D ex P 


/ 


iVoA 


PplP k ,D 


\ 




(27) 


and U = T 5 W with W denoting the transmission bandwidth. 

Proof: Let Y = maxji/j}^. Then, the complementary CDF of Y, F Y (x) = Pr[Y > 
x] — 1 — F Y (x), becomes fyM = 1 -n.^^.w. Hence, first recall the channel coherence 
consistency within the duration of sensing time. Then, according to ©, while using (fl8l) and 
(l22l) . (l26l) is derived in a closed formulation. ■ 

At this point it should be mentioned that a similar approach (but with full-blind relays instead) 
was studied in m, whereas assuming a single-sampled spectrum sensing, yielding a formula 
for the detection probability in terms of infinite series representation. 


D. Average Harvested Energy 

It is meaningful to investigate the average harvested energy of the ith relay (1 < i < M), 
given a fixed harvesting phase duration (c.f., Fig. [2]). 

Proposition 2: The average harvested power of the ith relay, defined as E H ,i, which is collected 
upon a detection of the primary system’s transmission, is given by 

Eh,% = Pa(X)E H ,u (28) 


where 


/ 


E H ,i = VPp f r 


r =1 k =1 


\ 


n 


while -Pd(A) is given by ([26 


7 P k ,Ri - lPj,Ri 

\rt k 


7p k ,Ri- 


(29) 


Proof: It holds that E = E [E H ,i] = ./7° xfE H , (x)dx, while using ([5]) and (fT71) . we have 


that 


L r 

1phJ X ) = fr 


r= 1 fc=l 




IT_ l _ 

7 j (lPk,Ri - 7 Pj,Ri) 




exp 


riPdPplP k ,Ri 




VPd(\)Pp 


Then, after some simple manipulations, (1281) is obtained. 


6 It is noteworthy that the corresponding average harvested power of the secondary source is directly obtained by substituting 
subscript Ri with S into d28t . denoting the corresponding node. 
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E. Outage Probability in the Transmission Phase 

Following similar lines of reasoning as in the reporting phase, while based on © and ([7]). the 
corresponding e2e SNR of the S — R s — D link is given by 


jr)jr> 

(T) _ 7l,s 12,s 
7e2 e,s ~ (T) 


7 ZJ+U 


(T) 


(30) 


where 


(T) A PS,R s |, 

7l, s = \ h S, 


C T ) 

2 , (T) a Pr s ,D 

Rs I and 72,s = 


h 


Rs,D 


(31) 

T , denoting 


N 0 1 ™ N 0 

Also, UP is explicitly defined in ([2Tb . by substituting the superscript (■)' R ' with 
the transmission phase this time. 

In what follows, it is important to investigate the scenario of asymmetric channels for the 
secondary nodes (i.n.i.d. statistics), in a sense that relays usually keep unequal distances between 
the two secondary end nodes (i.e., S and D) in practical dual-hop network setups. However, the 
scenario of symmetric channels (i.i.d. statistics), i.e., when the distances between the source-to- 
relay and relay-to-destination links are equal for each relay, is also considered, which can serve 
as a performance benchmark. 

1) Outage Probability: Outage probability, P out , is defined as the probability that the SNR of 
the e2e S — R s — I) link falls below a certain threshold value, 7 *, such that 


Pout (7th) = Pr 


(T) ^ 
7e2e,s < 7th 


(32) 


Lemma 2: CDF of the e2e SNR for the S — R s — D link over Rayleigh fading channels is 
expressed as 


F at) (x) = 1 - 2 E m exp 

'e2e,s 


N 0 x 


PS,RslS,R s J V ZmPs,Rs 1 S,R. 
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with 
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Proof: The proof is provided in Appendix |Bj 


(35) 


(36) 


Proposition 3: Outage probability of the e2e SNR for the secondary system, during its trans¬ 
mission phase, is presented for i.n.i.d. and i.i.d. channels, respectively, as 


M 


f« Li) (Ta,) = Z Pr [ ii ‘] F 7 ' T> W 

‘ <e2 e.i 


i= 1 


with 


and 


Pr[i2i] = 1 - 


j M—1 


(Aq + PRi,DlRi,DZM-l) Zm-1 


(37) 


(38) 


Z; d> (7,h) = F ,r,(7,h) 

Ie2e 


( 39 ) 
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where F en (•) is given by (1331) . while F (r)(-) is obtained by dropping index-/' from the latter 

Te2 e,i %2e 

expression, reflecting to identical statistics. 


Proof: The proof is given in Appendix 0 


At this point, it should be stated that the previously derived expressions regarding the detection 
and outage probabilities along with the average harvested power are in a closed form including 
finite sum series, whereas they are exact. Hence, they are much more computationally efficient 
than existing methods so far (e.g., the rather demanding Monte Carlo (MC) simulations and/or 
manifold numerical integrations). 

F. Avoiding Saturation of Fixed Gain Relaying 

When the received signal at the first hop is slightly faded or the link-distance between the 
involved nodes is quite small, the aforementioned fixed gain may go into saturation ll36l §V]. To 
avoid such a scenario, we introduce the modified fixed gain, suitable for practical applications, 
used for both the reporting and transmitting phases correspondingly. Starting from the former 
phase and recalling ©, it should hold for G -& ti that 



Gn,i Pp 


(40) 


where JCn represents a certain parameter that should not be exceeded during the AF process at 
the i\h relay. A numerical solution of JC-r is provided in Appendix iDl Thereby, substituting 
with F-n, the modified fixed gain is obtained, which results in clipping at the amplifier 

whenever (l40l) is violated. Furthermore, it is straightforward to show that the modified gain for 
the transmission phase becomes Gr,i\fFf with JC T denoting the corresponding constant value 
used for this phase. The solution of /C 7 - can be obtained by following similar lines of reasoning 
as for the derivation of the aforementioned JCr. 


IV. Energy Consumption 


Motivated by the general interests towards green communications in emerging and future 
wireless systems, we deploy the results of the previous sections to analyze the average energy 
required per node in the proposed CR cooperative system. Subsequently, a necessary condition 
for minimizing this quantity is presented. 
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A. Average Energy Consumption 

The average energy consumed at the zth relay in each frame (c.f., Fig. [2]) is given by 

Mil 4 E S ,T S U + E Ttl T R U + (1 - P rf (A))Pr [R,] E t ,T d - E H ,T D (41) 

where E Rl , E Rl and Ep A are the average energy consumed at the sensing, reporting and 
transmission phase, respectively. In addition, T s , T R and T D denote the duration of the sensing, 
reporting and transmission phase, correspondingly. Recall that T^W = U stands for the time- 
bandwidth product, i.e., the number of samples required at the sensing phase. Notice from 
m that a transmission event occurs with a probability 1 — P d { A), whereas Pr [ R,, denotes the 
probability that the zth relay is selected for transmission, as provided in (1381) for i.n.i.d. and 
(1C.4I) for i.i.d. setups, respectively. On the other hand, as previously stated, all the relays enter 
into the harvesting phase with a probability P d { A) during T D , when a primary transmission is 
sensed. Notice that P d { A) is already included within E R ,i, by referring back to (1281) . 

The sensing energy can be considered identical for all the secondary nodes and, thus, (l4TT) is 
simplified to 

Ptoil 4 E s T s U + E r .{P r U + ((1 - P rf (A))Pr [R,] P Tjl - P d (X)E H ^ T D . (42) 

Regarding E s , it holds that ll37l E s = Pr x , where P R:r is the circuit power used to capture 
the received signal(s) power. Moreover, for E R)i , we have that Il37j| - |]39l E R/l = [> ' R ' ] n + Pt x , 
where P Tx is the circuit power used for signal transmission and p IRr) is given by (flOl) . In 
general, both P Rx and P Rx are quite small, since each AF relay does not perform decoding, 
which is usually a more power-consuming operation lf37ll . Similarly, E T ,i is obtained as It38ll . 
[f39l E t , = ji/Pp) + Pt x , where p' Ri ] D is provided in (fl3l) . Finally, E Hl is presented in (l28l) and, 
hence, the average energy consumption is obtained in a closed-form. 

B. Energy Consumption Minimization 

A fundamental requirement for the energy minimization problem is the provision of an 
appropriate time duration for the transmission of a certain amount of data bits. On the other 
hand, this duration should also satisfy an appropriate energy harvesting level (when primary 
transmission is detected), which minimizes the energy consumption. Based on the proposed 
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mode of (hybrid) operation, while closely observing (1421 ) . the tradeoff between sensing time and 
transmission/harvesting time duration is not trivial. 

First, we assume that T R is fixed, since only the transmission of a probe message is involved 
into the reporting phase. Then, a maximization of Ts (i.e., a higher U ) would reflect to a higher 
detection probability by capturing the transmission of primary nodes, at the cost of minimizing 
Td- In turn, a shorter To would reduce the amount of transmitted data (in the case of transmission 
phase), whereas it would also reduce the amount of harvested energy (in the case of harvesting 
phase). Thereby, maximizing T s results to a higher energy consumption. On the other hand, 
minimizing Ts is also not a fruitful option, since the detection probability is reduced and, hence, 
a potential primary transmission may not be captured accurately. In such a case, a packet collision 
could occur during data transmission, while causing unexpected co-channel interference to the 
primary service. 

Based on (f27l) and (l29l) . (l42l) becomes 


S T S U + E r ,T r U + (A(A) TsW Pr [f?,] E Ti - E H , (l - A(A) TsW )) T D 



= E S T 2 S W + E Rii T R T s W + (A(Af sW Pr [R t ] E Tj * - E H>i (l - A(Af sW )) (T - T s ) 

(43) 


where T A T lotal — T R is the remaining duration of each frame, which is assumed as a fixed 
constant, whereas T total is the total frame duration. 


Then, the considered energy consumption minimization problem is formulated as 



s.t. V [i) > V * 


(44) 


where V* denotes a target amount of data bits that should be transmitted during the transmission 
phase and T>® stands for the amount of transmitted data of the /1h relay within a given 
transmission duration, which is defined as 


V w A (l - P d (A))Pr [. Ri} RT d = A(A ) TsW R(T - T s ) 
= A(A) TsW Pr [. R t ] W log 2 (l + 7 2i) (T - T s ) 


(45) 


with R A VVlog 2 (l + 7 e 2 e,j) denoting the transmission rate (in bps), while is explicitly 
defined back in (l30l) . 
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Keeping in mind that ln(A(A)) < 0 (since A (A) < 1), it readily follows from (1431) that 


d 2 E 


<0 

total 


d 2 T s 


= 2 E S W - A(A Y syv (E H:l + E t ,Pr [i2 i ])Wln(A(A))(2 - (T - T S )W ln(A(A))) > 0 

(46) 


revealing the convexity of the objective function. In addition, to extract the hidden convexity of 
the included constraint of (l44l) . the following transformation is applied 


A(Af» w Pr[i? i ] Wlogj (l + 7 A) (T-T s )> V 


A(A )~ T S w T>* 
2(3’-T i g)WPr[ii i ] 

S. , 


1 - 72, £ 


(47) 


A X>'(») 


Hence, the original optimization problem in (l44l) can be reformulated in the standard form as 


p(0 


total 


min E, 

T S > 0 

s.t. D' w < o. 


(48) 


Also, it holds that 


<9 2 P /(i) = J i n 2 / 2 a / A(A)~ r * w P* _ A(A)- r * w P*ln(A(A)) V 

d 2 T s “ | 1 J VPr fa] W(T - T 5 ) 2 Pr [i^] (T - T s ) J 


ln(2)A(A)" TsW 


2D* 


2D*ln(A(A)) D*Wln 2 (A(A)) A 


Pr [Ri] W(T - T s ) 3 Pr {R t } (T - T s ) 2 Pr [Jfc] (T - T s ) J 


> 0. 


(49) 

Thus, both the objective and constraint functions are convex with respect to Ts, which implies that 
the Karush-Kuhn-Tucker (KKT) conditions are necessary and sufficient for the optimal solution, 
whereas a unique minimum value exists fiOll . Motivated by the convexity of the problem, we 
introduce the following Lagrangian multiplier, termed /i, into the equation: 

£ = /i > 0 (50) 


and we set 


dC 

dTs 


= 0. 


(51) 
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Solving this system, it follows that 

A (A) _T S VV 'D* 

_ 2 (T-T S) wp r i« i jA(A) TsW Pr [i?,:] (T - T S ) 2 W 
/i_ X>*ln(2)(l - (T - T s )Wln(A(A))) 

x [ e h , + 2 E S T S W + E R}i T R W - A(A ) TsW (E H , + E T>i Pr [R t })( 1 - (T - T s )Wln(A(A)))) . 

(52) 


Since p > 0, the following necessary condition for optimally minimizing the average energy 
consumption yields as: 


E H ,i + E T l Vr [Ri] < 


A(A )- t s w (E h , + 2 E S T* S W + E r ,T r W) 
l-(T-r*)Wln(A(A)) 


(53) 


where Tg denotes the optimal sensing time. 

Interestingly, the instantaneous data rate does not affect the latter condition for the efficiency on 
the energy consumption since (1531) is independent of r i^ i . Obviously, the detection probability 
(obtained from (l26l) in a closed-form) plays a key role to the overall performance of the secondary 
system and to EE. Also, notice that for a sufficiently high detection probability (i.e., very low 
miss probability), A(A) —>■ 0 + , which implies that the right-hand side (RHS) of (1531) takes 
extremely high values. In turn, for A (A) —)■ 1~, RHS of (1531) maintains low values. Hence, (1531) 
is more likely to occur for an increased detection efficiency and vice versa. In other words, from 
a green communications viewpoint, the harvested energy is greatly enhanced as compared to the 
energy consumed for data transmission, as long as a robust and quite accurate detection scheme 
is preserved in the considered system. 

Overall, it is important to mention that in practical applications certain licensed spectrum bands 
allocated for primary users may be idle quite often (the so-called spectrum hole effect |[4T|| ). In 
this case, the secondary system would operate in consecutive time-frames, utilizing the previously 
proposed opportunistic relay scheduling scheme. Hence, it is highly possible that one or more 
secondary relays may reach a minimum energy level, defined as because the system enters 

more the transmission phase rather than the harvesting phase at this time period. The value of 
can be considered as an application-dependent fixed parameter and/or predetermined by 
the system manufacturer. As an illustrative example, consider the scenario when the system 
realizes a spectrum hole and the jth secondary relay has reached . Then, upon the next 

data transmission phase, the secondary receiver utili z es (1531) for the jth relay. If this condition 
is satisfied, the jth relay takes place in the following relay selection process. On the other 





21 


hand, if the latter condition is not satisfied, the secondary receiver excludes the jth relay for 
a potential selection of the data transmission phase. Thus, the analytical formulation of Pr [i2j] 
(with 1 < i < M and i ^ j) is provided back in (l38l) by substituting M with M — 1 to denote 
the absence of the jth relay during the selection process. 

V. Numerical Results 

In this section, numerical results are presented and cross-compared with MC simulations to 
assess our theoretical findings. Each MC run was conducted over 10 6 RV trials. There is a perfect 
match between all the analytical and the respective simulation results and, therefore, the accuracy 
of the presented analysis is verified. Henceforth, for notational simplicity and without loss of 
generality, we assume a common time correlation coefficient, defined as p. Moreover, 77 = 0.35, 
N 0 = — 131dBm, W = 1MHz, the system carrier frequency is 2.5GHz, while, as previously 
mentioned, the number of sensing samples is defined as U = T S W. Also, the path-loss exponent 
is assumed fixed as a = 4, corresponding to a classical macro-cell urban environment lf27l Table 
2.2]. In addition, Pr [9j = 1] = Pr [6j = 0] = 0.5 is assumed, without loss of generality, which 
is in agreement with other research works, e.g., If39l ll42l . All the included link distances are 
normalized with a reference distance equal to 1km. In addition, for clarity reasons, we assume 
that the distance between the Zth primary node and secondary source equals the distance between 
the Zth primary node and zth relay and the distance between the Zth primary node and destination, 
i.e., dp h s = dp h R i = dp h D — dp, V l, 1 . In what follows (owing to the non-identical statistics of 
the included nodes), we consider the following link-distance scenarios for the primary nodes; 
for L = 1 let dp 1 G M + (in km), while for L > 1 it is assumed that 

dp l+1 — dp t + 0.01 V/ G {1, L — lj. (54) 

Following similar lines of reasoning, when the link distances per hop for the secondary nodes 
are non-identical, it is assumed that 

ds,R i+ 1 — ds,Ri + 0.005, 

d Ri+1 , D = d^D + 0.005 Vi G {1, M — 1}. (55) 

Figure [3] illustrates the detection probability for various distances between primary and sec¬ 
ondary nodes. Its performance is worse for higher A threshold values and/or the existence of 
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d p (km) 

1 


Fig. 3. Pd{ A) vs. dp 1 . Recall ( 1541 . for the scenarios when L > 1. Also, M = 1, U = 200, dp,D = 0.1km, -2- = 2dB, and 


fewer primary nodes, as expected. This occurs due to the fact that when more primary nodes 
are placed in the vicinity of the secondary nodes, transmitting using a relatively high power, 
their active presence is more likely to be detected and vice versa. Scanning the open technical 
literature, a generally accepted target on the detection probability is found to be P£{ A) = 90% 
(e.g., see [[43]]). Then, as an illustrative example, when primary nodes keep a minimum distance 
of 0.4km from secondary relays and/or destination, the latter target is satisfied for the given 
system settings. 

In Fig. [4] outage performance of the secondary system is depicted for various P max values. It 
is worth noting that the diversity order is always one regardless of the number of relays, which 
is in agreement with [[4]. Moreover, when very rapidly varying fading channels are present 
(e.g., p = 0.1), adding more relays does not alter the coding gain either (in fact, there is quite 
a marginal performance difference, which can be considered as negligible). Hence, the overall 
outage performance, i.e., both the coding and diversity gains, cannot be enhanced by maximizing 
M in such environments. On the other hand, coding gain is improved by adding more relays 
into the system, when semi-constant channel fading conditions are realized (e.g., p = 0.9). 

Figures \5\ [6] and [7] are devoted to the average energy consumed by using the proposed 
opportunistic strategy, from a green communications perspective. To this end, P t otai is numerically 
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Fig. 4. Pout(7th) VS. Pmax/No- Also, L = 2, d s ,i{ = d R , D = 0.1km, d Pl}R = 0.3km. d PltD = 0.4km, = 3dB, 

= 6dB, and & = 30dB. 


evaluated using (l42l) under different distances dp, in Fig. [5] Obviously, there is an emphatic 
energy gain (on average) for each secondary node as this distance is relatively small (e.g., when 
d Pl < 0 ,4km). Also, EE is greatly enhanced for higher L values (i.e., more primary nodes). 
This is a reasonable outcome since the overall average harvested energy is increased in such 
a scenario. In general, it can be seen that the average harvested energy is higher than the 
corresponding energy that is consumed for sensing, reporting and (potential) transmitting. This 
beneficial phenomenon stops holding only for the case when primary nodes are rather far-distant 
(e.g., when d Pl > 1.1km). 

In Fig. [H E/total is numerically evaluated for various durations of the sensing time T s . Two 
different approaches are cross-compared; the proposed hybrid one (entitled as harvesting ap¬ 
proach) and a suboptimal approach which does not support harvesting (entitled as non-harvesting 
approach). The latter approach, supports either secondary data transmission or an idle state for the 
included nodes in the case of primary detection. Hence, the lack of energy harvesting manifests 
itself since the right-most depiction of Fig. [6] reveals a certain amount of consumed energy. On 
the other hand, the negative scale of E'totai at the left-most depiction describes an overall EE and 
an opportunistically power-saving system setup. 

An insightful observation obtained from Figs. [5] and [6] is the fact that the presence of more or 
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Fig. 5. .Etotai vs. various distances of the primary node-1 and the secondary system (i.e., dp x ). Also, i.i.d. statistics between 
the S — R and R — D links are considered with dp : p> = 0.1km and four relays (i.e., M = 4). Other parameters are fixed as: 
Pmsa = p P = 20dBm, A = Q = 17dBm, Pp x = lOdBm, Pp x = 9dBm, T tota i = 100msec, Tp = 1msec, and Tg = 20msec 
(thus Td = 79msec). 



Fig. 6. total vs. various durations of the sensing time Tg. Also, dp 1 = 0.4km. Other system parameters are the same as in 
Fig. 0 


less relays at the secondary system does not dramatically affect the average energy consumption. 
This occurs because when the system enters into the harvesting phase, upon the detection 
of primary transmission(s), all the included relays switch to the harvesting mode. Since the 
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Fig. 7. ECG vs. various durations of the sensing time Ts. Also, dp 1 = 0.5km, ds,^ = dR lt D = 0.2km, and M = 1. Other 
system parameters are the same as in Fig. [5] 


harvested energy is the dominant factor for EE and power savings, the derived results are quite 
straightforward. As previously noticed, performance gain from the introduction of more relays 
into the system appears only at the transmission stage of the secondary communication. 

In Fig. [7] the energy consumed for sensing, reporting and transmitting data (given the trans¬ 
mission phase) is evaluated with respect to the harvested energy (given the harvesting phase). 
To this end, the energy-consumption gain is introduced, namely, ECG, which is defined, using 


ECG A E sTs + Er'O'rTsW + (1 - P d ( A))Pr [Rj] E t ,T d 

P d (\)E Hji T D 

As can be seen from Fig. [7] the presence of two or more primary nodes provides an emphatic 
difference on power savings as compared to the rather overoptimistic scenario with only one 
primary node (i.e., when L = 1). This occurs not only due to the maximization of power collected 
during the harvesting phase, but also due to quite an enhanced detection probability for L > 2. 

Finally, Table U illustrates the optimal Ts (which can be obtained with the aid of several nu¬ 
merical optimization tooled)- given the optimization problem in (1481) with its involved constraints. 


7 For instance, a global optimization algorithm that can be employed for numerical estimation of Ts is the NMinimize function 
within the MATHEMATIC A™ software package. 






26 


TABLE I 

Optimal Sensing Time Tg (in sec) Obtained From ( f48li 


M 

L 

1 

2 

3 

4 

1 


0.0873 

0.0815 

0.0658 

0.000115 

2 


0.0873 

0.0815 

0.00372 

0.000115 

3 


0.0873 

0.0815 

0.00372 

0.000115 

4 


0.0873 

0.0815 

0.00370 

0.000106 


= 7dB, and {-SjF, j^} = 30dB, ds ,r 1 = d Rl ,D = 0.5km, d Pl , Rl = d PliD = 1km, T = 100msec, T R = 1msec, 
and R = 100Kbps. 


It can be readily seen that Ts is greatly minimized for higher L and M values, especially when 
M > 2 and L > 3. This is a natural outcome gained from the spatial diversity of cooperative 
sensing by many relay nodes and, thus, from the rather efficient detection probability within a 
minimum sensing time. 


VI. Conclusions 

A new green cognitive dual-hop relaying system with multiple AF relays was proposed and 
its performance is analytically investigated. The novelty of this work relied on the introduction 
of a hybrid mode of operation for the involved secondary nodes, which opportunistically switch 
between data transmission and energy harvesting, correspondingly, upon the detection of primary 
network activity. New closed-form expressions regarding important system performance metrics 
were obtained over Rayleigh channel fading conditions, such as the detection probability, outage 
probability, and average harvested energy of each secondary node. Capitalizing on the derived 
expressions, the energy conservation of secondary nodes was further investigated, by modeling 
and analyzing the overall energy consumption minimization problem. A necessary and sufficient 
optimality condition for the proposed opportunistic strategy was presented, while some useful 
engineering insights were manifested. 

A natural extension of the proposed scheme could be the analytical performance evaluation 
and/or optimization when the secondary nodes utilize multi-band spectrum sensing/spectrum 
operation, which represents a challenging topic for future research. 
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Appendix 

A. Derivation of 4221). ( [271 ), ( 1271) . and ( 1251 ) 

CDF of the e2e SNR for the ith secondary relay reads as l44l Eq. (15)] 


F% (®) = / E 1 w (z + -—- ) Lm(y)dy. 

% 2 J,i Jo 7l ^ V y ' ^ 


(A.l) 


Based on (l20l) . F^-r.) (x) can be directly obtained from (fl 8 l) . by substituting D with Ri. Addi¬ 


tionally, f (k)(x) stems as in © by substituting p t J and 7 ^ with Pp^ D and 7 Ri>D , respectively. 


72, i 


Therefore, utilizing [[241 Eq. (3.471.9)] into (IA.1I) . (l22l) can be easily extracted. 
Further, (1211) can be rewritten as 


u m = 


(z + l) f(K){x)dx 


7i ,i 


-1 


(A.2) 


while f (k)(-) is directly obtained from (fT71) by substituting subscript I) with R,. Hence, utilizing 


7i 


Eq. (2.3.4.2)] into (IA.2I) . (1231) is obtained. 

Regarding the derivation of (l24l) . referring back to (flQl) . we have that 


No exp — 


Nqx 


Pmax~t R-B 


L(K) (x) = < 
'2 


Pmax'YRi ,D 

Ar Tirr 1 ( N 0 E i<lR i } x \ 

jV 0 E[g fl .]exp^- ^ ^ J 
QlRi.D 




(A.3) 


1 E't'TR;] ^ p 


Q_ 

max 


Hence, it yields that 


F m{x) — 1 — (1 - F ' (k), p 

12,i \ '2,i Kmax 


0)) l 1 ~ 


E 1TrJ 


= 1 — exp 




lRi,D 


(A.4) 


By differentiating (IA.4I) . the corresponding (unconditional) PDF of 7 /, 7 y is formed as in (© with 


(P) 

the yielded transmission power p R . D defined in 

Finally, since E [q Ri ] = f 0 °° xf qR (x)dx, while based on (fill) and 01 Eq. (15)], it holds that 


L L L L ( IT k 

Z=1 k=0 n±=l rik = l ' 1 

7T1 ^ ^ 77./g ^ Af 

Thus, after some simple algebra, (l25l) arises. 


exp — 


E 


x 


lRi,Pi 7-7 lRi,P nt 


(A.5) 
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B. Derivation of fliJl) 

First, the transmission power for the S — R and R — D links are determined, accordingly, 
based on (fl2l) and (fl3l) . Similar to the derivation of (l24l) . ps,n s and p^ D are obtained in (l34l) 
and (1351) . respectively. Moreover, the S — R link follows a conditional (given the selected relay 
R s ) exponential distribution yielding 

F %]( x ) = 1 - exp f -——) • (B.l) 

7 ’ V PS,R a lS,R a J 

Regarding the second-hop of the transmission phase, recall that 7 ^ D and have correlated 
exponential distributions with a corresponding conditional PDF given by PTl Eq. (31)] 


fjn | *(T> (x\y) = 

’R S ,D I iR St D 


exp 


x+p 2 a y 


0--p^PR s ,DlR a ,D 


0 - pV)PR a ,DlR s ,D 


Zpsy/xy 


(1 - pV)Pr s ,d 1 r s ,d ) 


(B.2) 


From the theory of concomitants of ordered statistics Il46l § 6 . 8 ], it holds that 


fjr> 0) = / / (r) ,-(r) (x\y)f m ( y)dy. 

IR S ,D /r» iRx.D\ iRx.D iRx.D 


7 i2s, D I , D 


(B.3) 


Thus, since R s is selected based on maximizing the SNR of the relay-to-destination link, /. 


~en l 

7 Rs,D 


becomes 


Lcr) (y) = 'P M exp(-$ M y), (B.4) 

7 Rs,D 

by following similar lines of reasoning as for the derivation of (IA.5I) . Recall that + Ar and <I> M 
are defined in (1331) . Substituting tlB.2l) and (IB .41) into (IB.3I) . while utilizing |[47] Eq. (2.15.5.4)], 
it holds that 


/ AT) (x) = S M exp (~Z M x). (B.5) 

IR a ,D 

Further, using the integral in (IA.1I) . but plugging (IB.ll) and (IB.5I) . (1331) can be extracted. Finally, 

POO 

ul T) = 0 + 1 Y 1 f$p s {x)dx 
Jo 

with being the exponential PDF, yielding (l36l) . 
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C. Derivation of f li 71 ) and f \J9\ ) 

The unconditional outage probability of the received e2e SNR in (f37l) can be obtained by 
averaging (1331) over all the relay selection probabilities, yielding 

M 

-Pout** d} (7th) = ( 7 th), (C.l) 

' ^ 'e2 e,i 

i =1 

where Pr [/?,,] is defined in (HTI) . Suppose that Ii, is selected, which means that 7 Ri)D is the 
largest SNR among all R — D links. Let R v = rnaexcept the selected one (i.e., 
7 Ri,D > 7 Rv,d)- Then, Pr [R,] can be written as 

Pr [Ri] = Pr [max{ 7 Rl;D ,..., Jr m _ u d} < 1r v ,d\1r u d\ 

/»oo 

/ '^'lR v ,D flRi,D {pt)dx (C.2) 

Jo 

where F lRv D (-) is the CDF of the maximum of M — 1 i.n.i.d. exponential RVs. Based on (IB.5I) . 
but substituting M with M — 1, it holds that 

f , Jr d (x) = 1- exp(—Z M _ix). (C.3) 

Rv ’ D Z M _ X 

Hence, substituting (IC.3I) and © into (IC.2I) . (l37l) is directly obtained. 

In the special case of i.i.d. Rayleigh fading channels for the secondary system, we drop index-/' 
to ensure identical statistics per hop. To this end, the selection of each relay is equiprobable, 
yielding 

Pr[fli ] = T- (C.4) 

Therefore, (l37l) becomes (l39l) . thus, completing the proof. 


D. Numerical Solution of /Cr 

When the channel gain conditions of the first hop result to a violation of (|4(1 . the amplifier is 
clipped to /C-r. Consequently, similar to Il36l Eq. (19)], the modified fixed gain of the ith relay 
is modeled as 


/ 





(wf R, )"\ t!? <% 

-pf —, otherwise 

7^+1’ 


(D.l) 
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A Ki 14^^ 

where % = —-1 is a threshold value in terms of the received SNR at the /1h relay, given 

directly from (J40]), satisfying that G R Jp p J2f=i 0j\9Pj,Ri\ 2 + N o) = Gr- Then, averaging (ID.II ), 
the (unconditional) modified fixed gain yields as 
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where the second equality arises by utilizing ll24l Eq. (3.352.2)]. 

Finally, since = 1/(U lV] N {] ), fC- R can be numerically calculated by equating ( ID. 21 ) with 
the reciprocal of (l23l) . resulting to the following transcendental expression 
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satisfying a fixed gain that avoids saturation of the amplifier. 
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